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ABSTRACT 

This  study  examines  the  effectiveness  of  elastomer  coatings  for  specular, 
acoustic  backscatter  reduction  in  low  frequency,  underwater  appl ications. 
Four  representative  elastomers  were  selected:  Butyl  8252,  Neoprene  W, 
Polybutadiene,  and  Thiokol  RD.  Their  coating  reflection  loss  for  CW  and 
impulse  incident  acoustic  plane  waves  was  examined  using  techniques  and 
computer  codes  developed  based  on  a  geometric  acoustics  approach.  The  per¬ 
formance  of  these  elastomer  coatings  was  examined  in  the  frequency  range 
from  10  Hz  to  10  kHz,  the  pressure  range  0  to  500  psig,  and  the  temperature 
range  0°  to  30°C;  ranges  commonly  encountered  in  oceanic  applications.  The 
results  indicate  that  at  low  frequencies  (order  100  Hz),  the  available 
elastomers,  applied  as  homogeneous  coatings  up  to  20  cm  thick,  are  not  suf¬ 
ficient  to  obtain  reflection  reduction  greater  than  6  dB  for  the  specified 
conditions.  Greater  specular  backscatter  reduction  may  be  achieved  by 
using  layered,  inhomogeneous  coatings,  but  a  more  relevant  approach  for 
these  low  frequencies  may  involve  the  examination  of  structural  resonant 
response  to  the  incident  acoustic  field  since  the  body  size  and  acoustic 
wavelength  are  of  the  same  order  of  magnitude  (Gaunard,  1977).  The  struc¬ 
tural  resonant  mechanism  of  backscatter  is  probably  dominant  at  low  fre¬ 
quencies. 
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1.  INTRODUCTION 

1 .1  Project  Objectives 

There  is  currently  some  interest  and  concern  in  the  (J.S.  Navy  regarding 
low-frequency,  active,  acoustic  ASW  systems.  If  low  frequency  (~  100  Hz) 
acoustic  pulses  were  to  be  utilized,  such  systems  might  allow  long-range 
detection  (because  of  negligible  sea-water  absorption)  and  they  might  be 
difficult  to  counter  with  an  active  system  onboard  the  target.  The  poten¬ 
tial  for  passive  countermeasures ,  such  as  anechoic  coatings,  has  not  been 
completely  ruled  out  however.  In  particular  for  anechoic  coatings,  several 
recent  preliminary  calculations  (see  below)  have  indicated  a  surprising 
capability  and  it  was  the  purpose  of  this  investigation  to  refine  these 
calculations  and  extend  then  to  other  possible  coating  materials. 

The  technical  objective  of  this  study  is  to  assess  the  effectiveness  of  low 
frequency  anechoic  coatings  on  target  strength  reduction  and  make  recommen¬ 
dations  concerning  their  development  and  utility.  In  particular,  acoustic 
reflection  characteri sties  of  a  homogeneous  viscoelastic  layer  (an  elasto¬ 
mer  coating)  at  low  frequencies  (on  the  order  of  100  Hz)  in  water  are  exam¬ 
ined.  The  effects  of  viscoelastic  material  properties  and  coating  thick¬ 
ness  on  the  reflection  response  are  analyzed.  Selected  elastomers  are 
examined  for  their  effectiveness  in  this  application  over  the  range  of 
environmental  conditions  to  be  encountered  using  computer  codes  developed 
to  calculate  the  continuous,  plane-wave,  acoustic  reflection  coefficient 
and  the  impulse  reflection  response  of  a  homogeneous,  plane  coating  over- 
lying  a  rigid,  perfect  reflecting  structure. 

1.2  Background  and  Scope 
1.2.1  Background 

Interest  in  low  frequency,  underwater  acoustic  absorbers  was  sparked  by 
results  of  a  study  by  G.  V.  Borgiotti  (1981)  which  indicated  that  at 
100  Hz  reflection  reduction  on  the  order  of  10  dB  was  possible  with  coating 
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thicknesses  on  the  order  of  30  cm  for  the  elastomer  Thiokol  RD.*  R.  L. 
G ran  (1982)  corrected  florgiotti's  calculations  and  found  that  reflection 
reduction  of  this  level  could  be  achieved  with  coating  thicknesses  on  the 
order  of  5  cm  (Figure  1-1).  After  more  detailed  study  of  the  typical  pro¬ 
perties  of  viscoelastic  materials,  it  has  been  found  that  the  encouraging 
conclusions  of  these  two  previous  studies  were  the  result  of  an  underesti¬ 
mate  for  the  Poisson  ratio  of  Thiokol  RD  which  is  about  0.495  rather  than 
0.44  taken  by  the  two  studies.  Unfortunately,  the  more  realistic  value  of 
v  =  0.495  results  in  much  poorer  reflection  reduction  as  discussed  later  in 
this  report. 

THICKNESS  (cm.) 


Figure  l-l.  Calculated  reflection  reduction  for  Thiokol  RD  with  v  =  0.40 
and  0.44  (from  Gran,  1982) 


*  An  experimental  vibration  damping  elastomer  developed  in  the  1940‘s  and 
50's  which  is  no  longer  available. 
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1.2.2  Scope 

This  study  effort  consisted  of  four  basic  steps.  The  first  step  was  the 
development  of  the  theoretical  background  and  basic  algorithms  to  solve  the 
problem  of  plane  wave  reflection  reduction  from  a  plane  layer  overlying  a 
rigid  perfect  reflector.  The  second  step  involved  obtaining  material  pro¬ 
perty  data  for  available  elastomers,  running  the  algorithms  with  the  data, 
and  developing  the  final  presentation  format  for  the  output  data  results. 
The  third  step  was  to  analyze  these  computed  results  and  parameter  varia¬ 
tions  for  the  selected  coating  materials.  The  fourth  and  final  step  was 
the  examination  of  application  and  analysis  alternatives  for  underwater, 
low  frequency  target  strength  reduction. 

The  task  of  obtaining  the  "complete"  low  frequency,  dynamic  viscoelastic 
properties  of  available  elastomers  for  the  code  calculations  proved  to  be 
most  difficult  and  time  consuming.  Although  the  dynamic  properties  of  a 
wide  variety  of  elastomers  were  available,  the  information  for  each  indivi¬ 
dual  elastomer  was  often  incomplete.  As  a  result,  certain  properties  had 
to  be  estimated  based  on  "typical"  characteristics  and  values  for  elasto¬ 
mers.  This  required  a  more  detailed  study  of  viscoelasticity  and  an  exami¬ 
nation  of  the  sensitivity  of  reflection  reduction  to  these  estimates  by 
variation  of  these  parameters. 

1.3  Cone! usi ons 

After  computing  the  acoustic  reflection  reduction  capabilities  of  Thiokol 
RD  (with  reasonable  estimates  of  its  material  properties)  and  three  other 
representative  elastomers,  our  conclusion  is  that  for  typical  oceanic  con¬ 
ditions  (0°-30°C)  and  reasonable  coating  thicknesses  (1-20  cm),  the  elas¬ 
tomers  considered  do  not  provide  sufficient  reflection  reduction  in  the 
single  homogeneous  layer  configuration  for  low  frequency  (order  100  Hz) 
acoustics.  The  reflection  loss  was,  in  all  cases,  no  more  than  6  dB  for 
coating  thicknesses  less  than  20  cm.  The  acoustic  wavelength  at  100  Hz  in 
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the  elastomer  layer  is  about  two  orders  of  magnitude  larger  than  a  reason¬ 
able  coating  thickness.  Therefore,  unless  the  acoustic  wavelength  in  the 
layer  can  be  significantly  reduced  by  reducing  the  sound  speed  in  the 
layer,  the  effect  of  viscoelastic  coatings  having  reasonable  thicknesses  on 
specular  acoustic  reflection  should  be  negligible  at  low  frequencies. 

These  conclusions  do  suggest  that  for  effective  low  frequency,  underwater 
anechoic  coatings,  the  available  elastomers  might  be  modified  to  provide  a 
lower  coating  sound  speed.  One  technique  for  accomplishing  this  is  by 
introducing  inhomogeneities  such  as  air  voids,  which  convert  energy  from 
dilatational  modes  into  shear  modes  for  which  the  wave  speed  is  much 
slower.  Also,  using  such  materials  in  coatings  with  a  more  complex  geo¬ 
metry  (layering  and/or  selective  distribution  about  a  body),  significant 
target  strength  reduction  may  be  achieved.  Such  an  approach  should  include 
modal  analysis  which  encompasses  the  total  structure/coating  system  rather 
than  just  a  simple  geometric  acoustic  analysis  since  resonant  interaction 
of  the  structure  with  the  acoustic  field  may  be  just  as  dominant  a  mechan¬ 
ism  at  these  low  frequencies  where  the  acoustic  wavelength  is  on  the  order 
of  the  structure  size. 

1.4  Report  Structure 

The  theory  and  analysis  used  to  develop  the  algorithms  for  calculating  the 
continuous  wave  reflection  coefficient,  the  reflected  impulse,  and  the 
reflection  response  to  an  arbitrary  signal  are  discussed  in  Section  2.  An 
examination  of  the  relevant  viscoelastic  properties  with  respect  to  their 
frequency,  temperature,  and  pressure  dependence  and  their  effect  on  acous¬ 
tic  absorption  follows  in  Section  3.  A  parametric  study  of  four  represen¬ 
tative  elastomers  is  then  presented  in  Section  4  in  which  the  effects  of 
coating  thickness,  temperature,  density,  and  dynamic  moduli  are  discussed. 
Lastly,  observations,  alternatives,  application  constraints,  and  a  summary 
with  recommendations  for  further  analysis  are  given  in  Section  5. 
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2.  THEORY  AND  IMPLEMENTATION 
2.1  Assumptions  and  Conditions 

The  acoustic  reflection  algorithms  were  developed  using  the  theory  of  geo¬ 
metric  acoustics  applied  to  plane  acoustic  waves  incident  on  a  plane,  homo¬ 
geneous  layer  overlying  a  rigid,  perfect  reflecting  foundation  (Figure 
2-1). 


INCIDENT  PLANE  WAVE  FIELD 


COATING  T  H 

S/ /  /  /  /  '/  7~7~7  7  f 

RIGID  FOUNDATION 


Figure  2-1.  Coating  geometry  for  the  study 
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The  assumption  of  an  infinite  plane  layer  neglects  diffraction  effects 
around  a  finite  body,  but  since  the  acoustic  wavelength  and  body  size  are 
within  the  geometrical  scattering  region,  (kr  >  2)  (Figure  2-2),  and  a 
coating  has  negligible  effects  on  diffraction,  the  assumption  is  quite 
valid  for  this  study  of  anechoic  coatings. 


Figure  2-2.  Scattering  from  a  rigid  cylinder  with  radius  r  (k  =  acoustic 
wavenumber)  (adapted  from  Skudrzyk,  1971) 


The  plane  layer  represents  an  elastomer  coating  and  the  geometry  results  in 
the  maximum  specular  reflection  possible  since  the  structure  is  approxima¬ 
ted  by  a  perfect  reflector.  Since  only  specular,  monostatic  backscatter  is 
considered,  only  normally  incident  plane  acoustic  waves  are  of  importance. 
As  a  result,  only  dilatational  modes  of  wave  propagation  within  the  homo¬ 
geneous  plane  layer  will  be  considered  since  conversion  to  shear  modes 
should  be  negligible  at  normal  incidence. 


l 
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2.2  Continuous  Wave  (CM)  Reflection  Coefficient 


The  plane  wave  (pressure  amplitude)  reflection  coefficient  for  a  plane 
viscoelastic  layer  overlying  a  perfect  reflecting  foundation  is  given  by 
the  equation 


where 


R  = 


*12  +  6 


■121) 


1  +  R 


12 


Hf~  n2 
1  * 


phase  term 


°1  '  °2 
!12  =  0.  f 


boundary  reflection  coefficient 


—  cos  9 
21 


V5 


-  sin2  e 


'jJ  ,  0 


acoustic  frequency,  incidence  angle 


pr  ci 


fluid  density  and  sound  speed 


h  >  f>2  *  c2 


layer  thickness,  density,  and  sound  speed 


A  viscoelastic  absorbing  layer's  sound  speed  (di 1 atati onal  wave  speed  in 
this  case)  C2  is  complex,  resulting  in  a  complex  reflection  coefficient  R, 
which  accounts  for  absorption.  A  detailed  formulation  of  the  expression 
for  R  as  a  function  of  frequency,  incidence  angle,  layer  thickness,  and 
material  properties  is  given  in  Appendix  A. 
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2.3  Reflected  Impulse  Response 

The  reflected  impulse  response  is  essentially  the  time  series  representa¬ 
tion  of  the  broadband  reflection  response  of  the  coating.  It  is  a  compact 
representation  summarizing  the  coating  response  over  a  large  range  of  fre¬ 
quencies.  Initially,  the  technique  was  considered  to  calculate  the  coating 
response  to  a  very  short  (msec)  duration  pulse  like  that  produced  by  a 
water  gun  which  does  not  produce  a  CW  signal  (Figure  2-3).  For  long  range 
propagation  in  which  high  frequencies  are  absorbed  and  deeper  sources  are 
utilized,  the  water  gun  (a  cavitation  device)  cannot  be  used  and  signal 
pulses  are  of  much  longer  duration  (order  tens  of  msec).  In  these  long 
range  appl ications,  the  CW  (frequency)  representation  of  the  coating 
response  is  more  relevant  than  the  impulse  (time)  representation.  The 
technique  developed  to  obtain  the  impulse  response  from  the  CW  results  does 
provide  a  method  of  calculating  the  coating  time  response  to  an  arbitrary 
pulse  or  signal  whose  transmitted  spectrum  is  known. 


WATER  GUN  PRESSURE  WAVEFORM 


TIME  (m  sec) 

Figure  2-3.  Water  gun  waveform  (from  Buoyoucos,  st  a?., 1980) 
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The  plane  layer  reflected  response  as  the  ratio  of  the  reflected  pressure 
amplitude  Prefl  to  the  incident  amplitude  Pinc  for  an  impulse  (sudden  on- 
off  pulse  of  unit  amplitude)  is  given  by 


■  1  *U) 


JcOT 

e 


d  CO 


where 


oQ  frequency  limit  of  R(m)  data 

x  =  7 —  fx  sine  +  z  cose  -  c.t) 

C1  1 

R('jo)  CW  reflection  coefficient  as  function  of  u 


R( -to)  =  R  (oo)  analytic  continuation  for  negative  frequencies 


x 


t 


distance  along  and  distance 
from  layer-water  boundary 

time 


Note  that  the  finite  frequency  limit  ,o0  of  the  integral  is  an  approximation 
of  the  exact  impulse  response  (when  u0  =  ®)  since  information  on  R(w)  at 
frequencies  above  ioQ  is  not  included.  The  frequencies  above  u0  are  essen¬ 
tially  "filtered"  out  of  the  incident  and  reflected  impulse.  In  practice, 
the  integral  is  evaluated  by  use  of  a  discrete  Fourier  Transform  (i.e.,  rr7 
algorithm)  which  represents  a  further  approximation  by  replacing  the  inte¬ 
gration  with  a  discrete  summation.  A  detailed  formulation  of  the  reflected 
impulse  response  expression  is  provided  in  Appendix  A  with  a  comparison  of 
computation  results  with  an  exact  solution  showing  good  agreement. 
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2.4  Reflected  Response  to  an  Arbitrary  Signal 

The  plane  layer  reflected  response  to  an  arbitrary  waveform  or  a  pulse¬ 
shaped  acoustic  signal  can  be  computed  by 

+W°  i 

Is(t)  =  /  S(u)  R(«o)  eltJT  du 
"uo 

where  R(u>)  is  the  complex  CW  reflection  coefficient  and  S(u)  is  the  complex 
spectral  decomposition  (or  Fourier  Transform)  of  the  incident  signal.  The 
spectral  decomposition  of  an  arbitrary  transmitted  signal  can  be  obtained 
by  various  signal  processing  techniques  (i.e.,  FFT)  so  that  the  reflected 
response  of  a  plane  layer  to  such  an  acoustic  signal  can  be  predicted  with 
this  technique.  Essentially,  the  incident  signal  S(u))  is  "modulated"  by 
the  plane  layer  response  R(«)  resulting  in  the  reflected  signal  response  of 
the  coating  to  an  arbitrary  incident  signal.  The  inverse  problem  is  also 
of  interest  in  which  a  known  transmitted  incident  signal  S(w)  and  reflected 
response  Is(x)  are  measured  to  find  the  coating  CW  reflection  response 
R(w). 
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3.  V ISCDF.LASTIC  PROPERTIES 

3.1  Background  and  Assumptions 

Viscoelasticity  represents  a  class  of  material  behavior  between  that  of  a 
viscous  liquid  and  an  elastic  solid.  The  assumption  of  linear  viscoelas¬ 
ticity  is  a  good  approximation  of  the  behavior  of  most  elastomers  in  typi¬ 
cal  stress-strain  conditions.  Linear  viscoelasticity  allows  the  modeling 
of  stress-strain  behavior  by  systems  of  linear  viscous  elements  (dashpots) 
and  linear  elastic  elements  (springs)  in  series  and  parallel.  For  our 
application,  the  assumption  of  linear  viscoelasticity  applies,  and  the  fol¬ 
lowing  discussion  is  based  on  such  a  concept. 

3.2  Definitions  and  Relations 

The  single,  most  important  material  property  of  a  viscoelastic  layer  in 
determining  its  effect  on  the  reflection  coefficient  is  the  complex  dilata- 
tional  wave  speed  c  of  the  elastomer  material  making  up  the  coating  layer, 
assuming  negligible  conversion  into  shear  wave  modes  near  normal  incidence. 
The  complex  sound  speed  is  dependent  on  the  density  and  complex  moduli  of 
the  material.  This  complex  sound  speed  can  be  calculated  from  the  other 
material  properties  by  the  standard  relations 

y  i  A  p 

2  =  M  =  *  +  2G  =  _  3  b  =  2(l-v)G 

P  ”  P  ”  P  ~  (l-2v)p 

where 

p  is  the  density 

M  is  the  complex  longitudinal  modulus 

X.  is  the  Lame  constant 

G  is  the  complex  shear  modulus 

K  is  the  complex  bulk  modulus 

v  is  Poisson's  ratio 
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The  complex  moduli  quantify  the  relationships  for  the  magnitude  and  phase 
delay  between  stress  and  strain  and,  therefore,  account  for  the  speed  and 
attenutation  of  sound  in  the  material.  Often,  the  term  loss  factor  or  loss 
tangent  is  introduced  to  describe  the  relation  between  the  real  and  imagin¬ 
ary  parts  of  a  complex  modulus.  For  example,  the  loss  factor  for  the  shear 
modulus  is 

Gi 

=  g-  so  that  G  =  Gr  (1  +  i6g) 
r 

where  Gr  =  Re {G }  and  =  Im{G}.  Unless  t>1e  complex  longitudinal  modulus  M 
is  known,  at  least  two  other  dynamic  properties  (i.e.,  v  and  G,  or  K  and  G, 
etc.)  must  be  known  in  order  to  calculate  the  complex  sound  speed  c.  The 
complex  sound  speed  and  attenuation  are  functions  of  frequency  since  the 
complex  dynamic  moduli  of  elastomers  are  functions  of  the  frequency  of 
stress  loading. 

3.3  Frequency  Dependence 

Usually,  the  available  information  on  frequency  dependent  complex  moduli 
for  elastomers  is  in  the  form  of  '"he  real  part  and  the  loss  factor  for  the 
complex  shear  nodulus  G.  For  this  reason,  the  frequency  dependence  of  the 
dynamic  shear  modulus  and  its  effect  on  sound  absorption  at  different  fre¬ 
quencies  will  be  discussed.  The  shape  of  the  complex  dynamic  shear  modulus 
curve  of  viscoelastic  materials  as  a  function  of  frequency  is  quite  general 
and  can  be  represented  by  data  for  unvulcanized  natural  rubber  (Figure 
3-1). 


The  dynamic  behavior  of  the  elastomer  can  be  separated  into  four  different 
categories  corresponding  to  the  different  frequency  regions  shown  in  Figure 
3-1.  At  very  low  frequencies,  viscoelastic  materials  behave  more  like  vis¬ 
cous  liquids  in  which  deformation  can  continue  at  a  fixed  level  of  stress 
(i.e.,  creep),  and  viscous  dissipation  is  the  chief  form  of  energy  absorp¬ 
tion.  This  type  of  behavior  is  characteristic  of  the  "flow  region". 
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LOG  (REDUCED  FREQUENCY,  cps) 


Figure  3-1.  Dynamic  shear  modulus  data  (Gr  =  real  part,  Gi  =  imaginary 
part,  5  =  G^/G_  loss  factor)  for  unvulcanized  natural  rubber 
at  25°C  (from  Davey  and  Payne,  1964) 

In  the  "rubbery  elastic  region"  (at  medium  frequencies),  elastomers  behave 
in  the  rubber-like  manner  that  we  associate  with  these  materials.  In  both 
the  flow  and  rubbery  elastic  regions,  the  dynamic  modulus  remains  relative¬ 
ly  constant. 

At  higher  frequencies,  the  dynamic  behavior  of  elastomers  passes  through  a 
"transition  region"  from  a  rubber-like  behavior  to  a  glass-like  be.iavior. 
In  this  transition  region,  the  complex  modulus  increases  very  rapidly  with 
frequency,  and  the  loss  factor  reaches  a  maximum.  It  is  within  this  fre¬ 
quency  range  where  the  maximum  absorption  for  the  material  occurs.  At  even 
higher  frequencies  above  the  transition  region,  elastomers  behave  in  a 
glassy  manner  more  like  an  elastic  solid.  In  this  "glassy  region,"  the 
loss  factor  drops  to  very  low  values,  and  very  little  energy  is  absorbed  by 
the  material.  The  real  part  of  the  shear  modulus  Gr  is  nearly  constant  and 
approaches  an  asymptotic  value  at  infinite  frequency  while  the  imaginary 
part  G^  affecting  absorption  decreases. 
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3.4  Temperature  and  Pressure  Dependence 

For  the  range  of  operating  pressures  of  interest  in  underwater  coating 
applications  (less  than  500  psi),  the  sound  speed  of  elastomers  is  rela¬ 
tively  insensitive  to  changes  in  pressure  (Figure  3-2).  By  comparison,  for 
the  range  of  operating  temperatures  0°-30°C,  the  sound  speed  and  dynamic 
moduli  of  elastomers  change  considerably  more.  The  effect  of  temperature 
changes  on  the  dynamic  modulus  of  an  elastomer  is  to  shift  the  modulus  vs. 
frequency  curve  in  frequency  while  preserving  its  shape.  In  fact,  dynamic 
modulus  information  can  be  represented  by  a  master  curve  vs.  reduced  fre¬ 
quency  (aju)  at  a  reference  temperature  TQ  and  two  empirical  shift  con¬ 
stants  Ca  and  C(j.  The  amount  that  the  modulus  curve  is  shifted  in  frequen¬ 
cy  (ay)  to  represent  data  at  another  temperature  T  is  computed  from  the  ULF 
(Williams,  Landel ,  Ferry)  equation  (Ferry,  19BD). 


HYDROSTATIC  PRESSURE 

Figure  3-2.  Pres  ure  dependence  of  the  sound  speed  (cp)  for  Butyl  8252  at 
three  different  temperatures  (from  Capps,  Thompson  and  Weber, 

1  rtoi  \  r  * 
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where  the  temperatures  are  in  units  of  degrees  Kelvin.  Normally,  an 
increase  in  temperature  causes  a  shift  of  the  modulus  curve  to  higher  fre¬ 
quencies  while  a  decrease  in  temperature  shifts  the  curve  towards  lower 
frequencies  (Figure  3-3). 


Figure  3-3.  Temperature  and  frequency  dependence  of  the  dynamic  shear 
modulus  (Gr)  and  loss  factor  (6)  near  the  transition  region 
where  Tt  and  w*.  are  the  transition  temperature  and  frequency 
(from  Snowdon 1968) 


3.5  Classification  of  Elastomers 

Elastomers  can  be  classified  as  either  high  loss  or  low  loss  depending  on 
the  application,  operating  temperature  and  frequency  range.  Comparing  the 
complex  dynamic  modulus  information  to  classify  a  material  as  either  high 
loss  or  low  loss  requires  comparison  of  both  loss  factor  and  the  real  part 
of  the  dynamic  modulus  within  the  temperature  and  frequency  range  of  inter¬ 
est.  The  fact  that  one  elastomer  has  a  higher  loss  factor  than  another  is 
not  sufficient  to  classify  it  as  a  higher  loss  material.  In  many  cases,  an 
elastomer  with  a  higher  loss  factor  than  another  can  exhibit  less  damping 
and  energy  absorption. 
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For  two  elastomers  with  equal  loss  factors  at  a  temperature  and  frequency 
of  interest,  the  one  with  the  lower  value  of  the  real  part  of  the  dynamic 
modulus  will  exhibit  greater  absorption  and  damping.  In  general,  a  low 
value  of  the  real  part  of  the  dynamic  modulus  with  a  high  value  of  loss 
factor  results  in  greater  damping  and  energy  absorption  (a  high  loss  mater¬ 
ial).  The  complex  modulus  (both  real  part  and  loss  factor)  of  the  elasto¬ 
mers  must  be  compared  for  the  specific  range  of  operating  temperatures  and 
frequencies  since  the  behavior  of  some  elastomers  is  extremely  temperature 
and  frequency  sensitive. 

The  longitudinal  modulus  M  is  normally  one  to  two  orders  of  magnitude 
greater  than  the  shear  modulus  5  for  viscoelastic  (rubber-like)  materials. 
This  observation  is  exemplified  by  one's  experience  that  rubber-like 
materials  are  much  easier  to  bend  (shear  deformation)  than  to  stretch 
(longitudinal  ieformati on)  and  by  the  fact  that  shear  waves  exhibit  much 
slower  propagation  speeds  (c^  =  ^  G/p)  than  dilatational  waves  in  visco¬ 
elastic  mat..ria!  s.  For  the  temperature  and  frequency  ranges  of  interest  in 
this  study  (0°-30°C)  and  order  100  Hz  acoustics,  the  common  elastomers  have 
values  of  the  magnitude  of  the  dynamic  shear  modulus  on  the  order  of  ID7 
and  19d  dyn/cn 2  while  the  magn  tude  of  the  dynamic  longitudinal  modulus  is 
on  the  order  of  10d  to  10iJ  dyn/cm2.  In  these  operating  conditions,  the 
density  and  sound  speed  of  these  elastomers  are  nearly  the  same  as  water 
(p  =  1.0  gm/cm"*  and  c  =  1.5  x  10b  cm/s)  so  that  their  characteri Stic  acous¬ 
tic  impedances  (pc)  are  well  matched  for  application  as  underwater  anechoic 
coatings  (this  is  not  true  for  applications  in  air  in  which  the  impedance 
mismatch  is  a  major  problem). 
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4.  STUDY  OF  FOUR  REPRESENTATIVE  ELASTOMERS 

4. 1  Properties  and  Selection 

Based  on  the  high/low  loss  criteria  discussed  in  the  previous  section  and 
the  availability  of  dynamic  material  property  data,  four  specific  elasto¬ 
mers  were  selected  to  study  their  effectiveness  in  low-frequency,  under¬ 
water  anechoic  coating  applications  and  their  sensitivity  to  environmental 
and  design  changes  (i.e.,  temperature  and  thickness).  The  four  materials 

considered  herein  are  Butyl  B252,  Neoprene  W,  Polybutadiene,  and  Thiokol 
RD.  These  were  selected  based  on  their  low  values  of  the  real  part  of  the 
dynamic  shear  modulus  and  high  values  of  loss  factor  at  frequencies  on  the 
order  of  100  Hz  in  the  temperature  range  0°-30°C  compared  to  other  avail¬ 
able  elastomers. 

The  material  property  data  used  in  this  study  were:  density  (p),  real  part 
of  the  dynamic  shear  modulus  (Gr),  shear  loss  factor  (6^),  and  real  part  of 
the  sound  speed  (cr).  In  order  to  obtain  longitudinal  wave  propagation 
information  from  the  shear  data  provided  (Gr,  Sq) ,  it  was  necessary  to 
assume  a  constant,  real  °oisson's  ratio  v  over  the  temperature  and  frequen¬ 
cy  range  of  interest.  This  assumption  probably  predicts  more  longitudinal 
absorption  than  is  actually  realized  since  it  results  in  the  longitudinal 
and  shear  loss  factors  being  equal  (5^  =  6q) .  At  the  1  iw  frequencies  of 

interest  here  (order  10D  Hz)  the  error  introduced  by  tfr  *  i  sumpti  -  hould 
not  affect  our  conclusions  appreciably.  Poisson's  ratios  for  tne  elasto¬ 

mers  were  calculated  from  cr,  Gr,  and  p  data  using  the  relation 

.  2( 1- v)G 

Cr  =  (l-2v)p 

assuming  that  c^  >>  c?  (which  is  reasonable  at  these  frequencies).  The 
results  of  these  calculations  are  shown  in  Table  4-1.  For  elastomers, 

Poisson's  ratio  is  typically  in  the  range  from  0.490  to  0.499.  A  reason¬ 
able  estimate  for  Thiokol  RO  is  then  v  =  0.495. 
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Table  4-1.  Density  (p),  sound  speed  (cr),  and  Poisson's  ratio  ( v) 

at  T  =  10°C,  f  »  102  Hz  for  the  four  elastomers  examined 


p(  gin/cm3) 

cr(x  105  cm/s) 

V 

Butyl  B252 

1.17 

1.58 

0.497 

Neoprene  W 

1.44 

1.60 

0.499 

Pol ybutadi ene 

1.13 

1.40 

0.498 

Thiokol  RD 

1.25 

1.13  * 

0.495  * 

*  unknown,  reasonable  estimate 


i 


The  master  curves  and  WLF  shift  constants  (C-,  and  Ch)  for  the  dynamic  shear 
modulus  were  available  for  Butyl  8252,  Neoprene  W,  and  Polybutadiene 
(Roberts  and  Madigosky,  1080;  and  Capps,  Weber  and  Thompson,  1981),  but  not 
for  Thiokol  00.  Data  for  Thiokol  RD  were  available  at  three  different  tem¬ 
peratures  (Snowdon,  1968)  and  were  interpolated  for  10°C.  The  different 
temperature  curves  then  allowed  calculation  of  the  WLF  constants  for  the 
master  curve  of  Thiokol  RD  at  10°C.  The  dynamic  shear  modulus  data  for  all 
four  materials  are  shown  in  Figure  4-1  followed  by  their  shear  loss  factors 
in  Figure  4-2. 


The  dynamic  longitudinal  modulus  data  calculated  from  the  results  in  Table 
4-1  and  Figure  4-1  using  the  relation 


are  presented  in  Figure  4-3.  The  longitudinal  wave  speed  used  in  the 
reflection  algorithms  can  then  be  found  from  the  relation 


« 
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FREQUENCY  (HZ)  FREQUENCY  (HZ) 


OHins 


PGLYBUT  RD 


(2»«W3/NjU3) 


t2»«H3/NJlO) 


Calculated  longitudinal  modulus  data  (at  10°C)  assuming  constant,  real 
Poisson's  ratio  v  and  6M  =  for  the  four  selected  elastomers  used  in 
computing  reflection  loss  ( -  real  part  Mf , - imaginary  part  Mj ) 


t-°l  ,-0I  ,.01  ,.01  ,.0l  ,.01  „ot 
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where  both  c  and  M  are  complex  functions  of  frequency  and  temperature. 
Another  expression  relating  shear  properties  to  the  longitudinal  properties 
is  given  by 

5M  =  2  5G  TT 

(Madigosky  and  Fiorito,  1982)  and  results  in  substanti all y  di fferent  values 
for  the  longitudinal  moduli  than  the  ones  calculated  by  assuming  a  con¬ 
stant ,  real  Poisson's  ratio.  This  expression  results  from  ultrasonic  data 
on  liquids  and  polymers  and  may  not  apply  at  the  low  frequencies  in  our 
case.  Essentially  this  high-frequency  expression  results  in  negligible 
longitudinal  absorption  (reflection  loss  less  than  1  dB)  since  Mr  is  about 
two  orders  of  magnitude  larger  than  Gr  for  elastomers.  The  longitudinal 
moduli  were  calculated  using  this  ul trasoni c-extrapol ated  expression  and 

the  results  are  displayed  in  Figure  4-4  followed  by  the  calculated  6j.  in 

Figure  4-8.  The  nearly  constant  value  of  Mp  over  frequency  (Figure  4-4) 
leads  us  to  believe  that  this  ultrasonic  extrapolation  is  not  applicable  to 

the  present  low  frequency  case  since  Mr  should  increase  with  frequency  in 

this  region.  Therefore,  the  longitudinal  data  calculated  using  the  con¬ 
stant,  real  Poisson's  ratio  assumptions  (Figure  4-3)  is  used  to  compute 
reflection  loss  in  this  study  and  should  be  valid  in  the  frequency  range  of 
i nterest. 

4.2  Parameters  and  Conditions 


For  low  frequency,  underwater,  anechoic  coating  applications  the  environ¬ 
mental  and  design  parameters  of  interest  are: 

acoust’c  frequency 

temperature 

pressure 

coating  thickness. 
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Typical  ranges  of  values  for  these  parameters  in  this  application  are  given 
in  Table  4-2.  The  ambient  water  density  and  sound  speed  were  taken  as  1.0 
gm/cmJ  and  1500  m/s,  respectively,  and  the  ambient  pressure  was  set  at  0 
psig  since  the  elastomers  were  found  to  be  quite  insensitive  to  pressure  in 
the  range  of  interest  (Figure  3-2).  The  effects  of  temperature  in  the 
range  0°-30°C  and  coating  thickness  in  the  range  1-20  cm  on  effective  coat¬ 
ing  reflection  reduction  at  100  Hz  were  examined  in  this  study.  The  dyna¬ 
mic  viscoelastic  properties  of  the  four  elastomers  examined  at  100  Hz  are 
tabulated  in  Table  4-3. 


Table  4-?.  Parameter  value  ranges  of  interest  for  underwater 
anechoic  coating  applications 


acoustic  frequency 

order  100  Hz 

temperature 

0°  -  30°C 

pressure 

0  -  500  psig 

coating  thickness 

1-20  cm 

Table  4-3.  Dynamic  modulus  data  for  the  four  selected  elastomers 

at  100  Hz  for  10°C  and  20°C. 


Gr  (10; 

dyn/cm2) 

Mr  (109 

dyn/cm2) 

6g  or 

10°C 

20°C 

10°C 

20°C 

10°C 

20°C 

Butyl  B252 

5.7 

2.9 

9.5 

4.9 

0.96 

0.88 

Neoprene  W 

3.1 

2.4 

15.0 

12.0 

0.25 

0.15 

Pol  yb'itadiene 

2.3 

0.9 

5.8 

2.3 

0.70 

0.23 

T'niokol  R0 

12.0 

2.2 

12.0 

2.2 

1.30 

1.21 
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4.3  Reflection  Loss  Results 
The  reflection  loss  is  defined  as 

RL  -  -10  log  | R ) 2  (dB) 

where  R  is  the  complex  CW  pressure  amplitude  reflection  coefficient  dis¬ 
cussed  in  Section  2.2.  The  lower  limit  of  RL  is  0  dB  corresponding  to 
total  reflection  of  incident  acoustic  energy.  The  larger  RL  (in  dB) ,  she 
less  acoustic  energy  is  reflected  and  the  greater  the  absorption  by  the 
coating.  It  is  well  to  reiterate  that  the  results  presented  are  for 
specular  backscatter  from  a  plane,  viscoelastic  coating  layer  on  a  perfect 
reflecting,  rigid  half-space  in  water. 

At  10°C,  all  four  selected  elastoners  exhibit  little  reflection  loss  at 
frequencies  below  1  kHz  (Figure  4-6).  Neoprene  W  and  Thiokol  RD  exhibit 
the  lowest  calculated  reflection  losses  at  this  temperature  due  to  their 
high  values  of  Mp  (in  spite  of  the  relatively  high  value  of  for  Thiokol 
RD;  see  Table  4-3).  Polybutadiene  exhibits  the  most  reflection  loss  in 
this  case  because  if  has  the  lowest  Mr  of  the  selected  materials  at  this 
temperature  (Table  4-3). 

At  2D°C,  the  CW  reflection  loss  results  are  quite  different  (Figure  4-7). 
Thiokol  RD  now  exhibits  the  most  reflection  reduction  since  it  has  the  low¬ 
est  Vr  and  highest  5m  at  this  temperature  (Table  4-3).  The  results  for  the 
other  three  elastomers  are  not  significantly  different  from  those  at  10°C 
except  that  Butyl  B252  now  exhibits  slightly  more  reflection  loss  than 
Polybutadiene  in  this  case.  At  all  temperatures  of  interest,  Neoprene  W 
gives  the  lowest  reflection  loss  of  the  selected  materials  because  it  has  a 
significantly  higher  Mp  and  lower  6^. 
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Comparing  the  CW  reflection  loss  results  at  the  two  temperatures  (10°C  and 
20JC) ,  it  is  apparent  that  the  selected  materials  exhibit  much  greater 
reflection  loss  at  high  frequencies  above  1  kHz  (except  for  Thiokol  RD  at 
iO^C).  In  fact,  maximum  reflection  reduction  for  these  elastomers  occurs 
in  the  neighborhood  of  If)  kHz  corresponding  to  the  quarter-wavel engtn 
thickness  frequency  in  which  the  elastomer  coating  has  a  "optimum"  thick¬ 
ness  to  maximize  reflection  loss.  41so,  Thiokol  RO  is  extremely  tempera¬ 
ture  sensitive,  changing  from  a  low  loss  elastomer  at  10°C  to  a  high  loss 
elastomer  at  20°C.  The  other  three  elastomers  were  relatively  insensitive 
to  temperature.  Figure  4-3  displays  the  temperature  and  thickness  depen¬ 
dence  of  coating  reflection  loss  at  100  Hz  for  the  four  elastomers.  It 
appears  that  Thiokol  RD  has  an  optimum  coating  thickness  near  13  cm  at  30JC 
while,  at  lower  temperatures ,  the  reflection  loss  increases  with  coating 
thickness.  For  Polybutadiene  at  100  Hz,  the  optimum  temperature  is  10°C 
since  the  reflection  reduction  is  less  at  both  lower  and  higher  tempera¬ 
tures  (Figure  4-3). 

Classifying  the  selected  elastomers  according  to  reflection  loss  at  100  Hz, 
Thiokol  RO  is  a  hign  loss  elastomer  at  temperatures  greater  than  20°C  and  a 
low  loss  elastoner  for  temperatures  below  10°C.  Butyl  32S2  is  a  medium 
loss  elastoner  at  20°C  and  above  and  a  low  loss  elastomer  at  10°C  and 
below.  Both  Polybutadiene  and  Neoprene  W  have  very  poor  reflection  reduc¬ 
tion  characteristics  and  can  be  classified  as  low  loss  elastomers.  None  of 
the  elastomers  provided  more  than  about  6  dB  of  reflection  reduction  Tor 
coating  thicknesses  less  than  20  cm  (Figure  4-8). 

The  sensitivity  of  the  calculated  reflection  loss  to  errors  in  the  estima¬ 
ted  Poisson's  ratio  v  is  examined  in  Figure  4-9.  Remember  that  Poisson’s 
ratio  for  elastomers  in  taese  conditions  is  within  the  range  from  0.490  to 
0.499.  4s  shown,  reflection  reduction  is  quite  sensitive  to  Poisson  s 
ratio.  The  sensitivity  of  reflection  reduction  to  density  was  found  to  be 
negligible  over  the  range  from  1  gm/cm3  to  2  gm/cm3. 
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4.4  Impulse  Reflection  Results 

The  reflection  of  an  imouVse  (unit  on-off  pulse)  from  plane  layers  of  the 
four  selected  elastomers  at  10°C  and  20°C  was  calculated  using  the  complex 
pressure  reflection  coefficient  vs.  frequency  results  and  the  technique 
discussed  in  Section  2.3.  The  CW  reflection  coefficient  (R)  data  covered 
frequencies  up  to  10  kHz,  so  for  our  calculations,  the  incident  impulse 
(Figure  1-10)  does  not  have  energy  at  frequencies  above  10  kHz.  The 
impulse  reflection  results  are  then  expressed  as  a  pressure  ratio  (the 
ratio  of  reflected  pressure  to  incident  impulse  pressure)  vs.  time  (Figures 
4-11  and  4-12)  or  time  series  of  the  reflected  pulse. 


-.004  -.002  0  .  002  .  004  .006  .008  .010 

7 DC  (SEC) 


Figure  4-10.  Incident  unit  impulse 
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Figure  4-11.  Calculated  impulse  reflection  from  a  1')  cm  layer  of  the  selected  elastomers 
at  10°C.  The  peak  for  Thiokol  RD  is  wide  due  to  FFT  resolution  (only 
components  up  to  1  kHz  are  included) 
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Since  the  elastomers  exhibit  larger  reflection  reduction  at  frequencies 
above  1  kHz,  the  impulse  reflection  results  are  much  better  than  for  100  Hz 
CW.  The  great  difference  between  the  results  at  10°C  and  20°C  for  Thiokol 
RD  clearly  demonstrates  the  large  temperature  sensitivity  of  its  reflection 
response.  The  other  three  elastomers  do  not  exhibit  nearly  as  much  temper¬ 
ature  dependence  in  the  impulse  reflection  results.  Note  that  these 
results  are  for  an  ideal  incident  impulse  (with  frequency  components  above 
10  kHz  filtered  out)  rather  than  a  realistic  impulse  which  has  some  charac¬ 
teristic  frequency  and  time  width.  The  reflection  of  a  realistic  impulse 
or  signal  of  arbitrary  shape  can  be  calculated  using  this  technique  (dis¬ 
cussed  in  Section  2.4).  These  impulse  reflection  results  provide  a  measure 
of  the  total  broadband  reflection  response  of  a  coating,  and  they  exemplify 
the  capaoility  to  calculate  the  coating  time  response  to  an  arbitrary 
acoustic  signal  using  this  technique. 

4.5  Limitations  of  the  Study 

The  major  limitation  on  the  present  study  was  the  availability  of  elastomer 
dynamic  properties.  This  limitation  was  overcome  by  combininq  available 
data  to  obtain  reasonable  estimates  for  the  unknown  material  properties 
(i.e.,  Poisson's  ratio)  at  the  temperatures  and  frequencies  of  interest. 
Four  elastomers  were  chosen  to  provide  a  cross-section  of  commonly  avail¬ 
able  materials  that  might  be  used  as  anechoic  coatings.  Both  Neoprene  W 
and  n0i ybutadi ene  represent  low  loss  elastomers  in  our  conditions.  At  tem¬ 
peratures  above  20°C,  Butyl  B252  represents  a  medium  loss  elastomer,  and 
Thiokol  RD  represents  a  high  loss  elastomer.  At  temperatures  below  10°C 
all  four  elastomers  can  be  considered  low  loss. 

The  dynamic  material  properties  available  were  limited  to  frequencies  up  to 
10  kHz  (for  the  temperature  range  of  interest)  so  the  effects  of  absorption 
at  higher  frequencies  are  not  considered.  This  is  a  good  apnroximation  for 
long  rang;  propagation  in  the  ocean  where  high  frequency  acoustic  wave  com¬ 
ponents  are  attenuated  significantly.  Also,  in  oceanic  applications  at  low 
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frequencies  on  the  order  of  100  Hz,  the  complete  backscatter  response  of  a 
coating  on  a  structure  whose  size  is  on  the  order  of  an  acoustical  wave¬ 
length  may  not  be  accurately  modeled  by  the  response  of  an  infinite  plane 
layer  because  other  backscatter  mechanisms  may  be  important.  Instead,  this 
study  examines  the  reflection  reduction  effectiveness  of  various  types  of 
elastomers  in  the  specular  reflection  backscatter  regime.  The  plane  layer 
reflection  is  a  good  measure  of  the  specular  reflection  reduction  of  these 
elastomers  for  application  in  anechoic  coatings  regardless  of  structure 
size  or  shape.  The  study,  therefore,  is  most  useful  for  applications  such 
as  material  selection  where  specular  reflection  reduction  by  a  homogeneous 
coating  is  of  primary  concern. 
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5.  DISCUSSION  AND  RECOMMEND ATIONS 

5.1  Alternatives  and  Options 

From  the  results  of  this  study,  we  find  that  elastomer  coatings  used  in  the 
simple  homogeneous ,  single  layer  confi gurati on  do  not  provide  significant 
low  frequency  acoustic,  specular  backscatter  reduction  for  coating  thick¬ 
nesses  of  20  cm  or  less  in  the  ocean.  At  low  frequencies  (on  the  order  of 
100  Hz)  in  the  ocean,  the  acoustic  wavelength  is  two  orders  of  magnitude 
larger  than  the  typical  coating  thickness.  Since  the  longitudinal  wave 
speeds  in  the  elastomers  are  almost  the  same  as  the  sound  speed  in  the 

water- ,  the  elastomer  coatings  are  incapaole  of  providing  significant 
absorption  at  these  wavelengths.  In  order  to  increase  the  absorption  of 
incident  acoustic  energy,  the  acoustic  wavelength  within  the  coating  can  be 
reduced  by  mode  conversion  of  longitudinal  wave  energy  into  shear  wave 
energy  for  which  the  propagation  speed  is  much  slower.  This  effectively 
reduces  the  coating  sound  speed  so  the  acoustic  energy  effectively  "sees" 
more  coating  and  absorption  as  it  propagates  at  the  slower  speed. 

One  way  to  achieve  a  lower  sound  speed  in  the  elastomer  coating  is  to 

introduce  inhomogeneities  such  3s  air  voids  which  convert  acoustic  energy 
to  shear  nodes  of  propagation.  There  is  a  trade-off  involved  in  this 
approach.  As  the  effective  sound  speed  in  the  coating  is  reduced,  the 

coating  absorption  capability  is  increased,  but  the  impedance  (pc)  mismatch 
between  the  water  and  coating  is  increased.  This  impedance  mismatch  causes 
a  portion  of  the  acoustic  energy  to  be  reflected  directly  from  the  water¬ 
coating  interface  without  propagating  (and  being  absorbed)  in  the  coating. 
Further,  such  coating  designs  may  be  more  sensitive  to  pressure  changes. 
Clearly,  there  is  some  optimum  coating  density  and  sound  speed  that  mini¬ 
mizes  the  total  acoustic  energy  reflected  from  both  the  water-coating 
boundary  and  the  propagation  through  the  coating  and  reflection  from  the 

structure  (Scharnhorst ,  1979) . 
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Multi-layer  coatings  may  also  help  to  reduce  the  impedance  mismatch  problem 
of  introducing  voids  by  matching  the  pc  impedance  at  the  water-coating 
boundary  and  allowing  a  transition  to  lower  sound  speed  and  greater  absorp¬ 
tion.  Such  an  approach  has  been  examined  at  high  frequencies  (greater  than 
1  kHz)  and  found  to  be  quite  effective  (Madigosky,  1981),  but  its  effec¬ 
tiveness  at  low  frequencies  (below  1  kHz)  is  still  in  question.  The  high 
frequency,  multi-layered,  inhomogeneous  elastomer  coating  approach  relies 
chiefly  on  quarter  wavelength  attenuation  within  the  layers  to  maximize 
reflection  reduction.  Comparable  levels  of  anechoic  effectiveness  should 
not  be  expected  at  low  frequencies  since  the  layer  sound  speed  would  have 
to  be  reduced  to  about  100  m/sec  or  less  for  a  20  cm  thick  layer  in  order 
to  take  advantage  of  the  quarter  wavelength  effect  at  100  Hz. 

Instead,  at  low  frequencies,  an  analysis  of  the  structure-coating  system 
interaction  with  the  acoustic  field  in  the  water  may  be  more  appropriate 
since  the  acoustic  wavelength  in  water  is  on  the  order  of  the  total  struc¬ 
ture  size  and  specular  reflection  might  not  be  the  primary  backscatteri ng 
mechanism.  In  fact,  resonant  interaction  of  the  coating-structure  and  the 
acoustic  field  may  cause  backscatter  greater  than  that  predicted  by  specu¬ 
lar  reflection  alone  (Gaunard,  1977).  There  are  indications  that  at  these 
acoustic  wavelengths,  the  dominant  mechanism  for  backscatter  is  not  specu¬ 
lar  reflection  but  resonant  structural  vibration  and  reradiation  by  the 
structure  (Gaunard  and  Ka.nir.s,  1982).  If  this  is  the  case,  a  modal  analy¬ 
sis  approach  to  examining  the  effects  of  low  frequency  coatings  is  more 
appropriate  than  the  geometric  acoustics  approach  of  this  study.  The  coat¬ 
ing  is  then  no  longer  an  "anechoic"  coating  in  the  sense  of  reducing 
reflection,  hut  serves  more  to  modify  the  modal  response  of  the  structural 
vibration  and  reradiation. 
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5.2  Constraints  for  IJndarwater  Applications 

In  the  pressure  range  for  underwater  anechoic  coating  applications  (0-500 
psig),  the  common  homogeneous  elastomers  are  essentially  insensitive  to 
pressure  changes.  Instead,  the  primary  sensitivity  to  depth  changes  is  due 
to  changes  in  seawater  temperature  (within  the  range  0°-30°C).  The  elasto¬ 
mers  with  broad  transition  regions  and  broad  loss  factor  vs.  frequency 
peaks  tend  to  be  less  temperature  sensitive  than  those  with  sharper  loss 
factor  vs.  frequency  peaks  (i.e.,  Thiokol  RD).  This  means  that  the  elasto¬ 
mers  which  tend  to  be  less  temperature  sensitive  also  tend  to  be  less  fre¬ 
quency  sensitive  and  would  be  advantageous  for  undersea  applications  if 
they  can  provide  significant  reflection  reduction.  Unfortunatel y ,  the 
elastomers  examined  do  not  provide  significant  reflection  loss  at  frequen¬ 
cies  below  1  kHz  throughout  the  range  of  application  conditions. 

If  the  coating  thickness  were  not  limited,  sufficient  anechoic  response 
could  be  accomplished  but  only  at  much  greater  thicknesses  than  the  range 
examined  (1-20  cm).  The  permissible  coating  thickness  is  limited  mainly  by 
the  coating  density  (buoyancy  considerations)  and  hydrodynamic  considera¬ 
tions.  'dost  common  elastomers  have  densities  slightly  greater  than  water, 
and  if  voids  were  introduced,  their  densities  would  decrease.  If  buoyancy 
was  the  only  factor  limiting  coating  thickness,  very  thick  inhomogeneous 
(void)  coatings  would  be  permissible.  The  effects  of  the  added  bulk,  mass, 
and  compliance  on  hydrodynamic  performance  would  have  to  be  examined  for 
such  thick  coatings.  It  appears  doubtful  that  any  available  elastomer  used 
as  a  homogeneous  coating  of  the  thickness  considered  would  provide  suffi¬ 
cient  specular  reflection  reduction  at  low  frequencies  over  the  temperature 
range  of  interest.  Thiokol  RD  does  provide  adequate  reflection  loss  at 
20°C  and  above,  but  it  is  extremely  temperature  sensitive  and,  therefore, 
not  suitable  for  these  applications  unless  the  coating  temperature  can  be 
maintained  above  20°C  (by  utilizing  waste  heat  from  the  vehicle's  power 
plant,  for  example).  If  maintaining  the  coating  temperature  is  considered, 
availability  of  sufficient  thermal  flux  would  have  to  be  examined. 
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5.3  Summary  and  Recommendations 

In  summary,  we  have  found  that  the  commonly  available  elastomers  are  quite 
ineffective  in  reducing  specular  reflection  at  low  frequencies  (order  100 
Hz)  for  the  tenperature  range  of  interest  (0°-30°C).  Thiokol  RD  does  pro¬ 
vide  adequate  reflection  loss  if  the  coating  temperature  can  be  maintained 
above  20JC,  perhaps  by  utilizing  waste  heat  (from  exhausted  power  plant 
cooling  water).  Introduction  of  more  complex  coating  geometries,  such  as 
inhomogenei ties  and  layering,  may  improve  coating  effectiveness  in  specular 
reflection  reduction,  but  if  the  acoustic  wavelength  is  on  the  order  of  the 
body  size,  resonant  response  of  the  whole  structure  system  may  be  the  domi¬ 
nant  backscatter  mechanism. 

The  geometric  acoustics  analysis  used  in  this  study  for  specular  reflection 
reduction  is  not  relevant  to  the  case  where  structure  vibrational  response 
is  the  dominant  backscatteri ng  mechanism.  Instead,  a  modal  analysis 
approach  to  the  total  structure  vibrational  response  with  consideration  of 
coating,  structure,  and  internal  fluid  properties,  geometries,  and  coupling 
would  be  necessary.  The  rinjor  mechanism  in  reducing  backscatter  in  this 
case  would  be  vibration  damping  rather  than  coating  absorption  of  incident 
acoustic  energy  as  in  the  case  of  specular  reflection  reduction  examined  by 
this  study. 
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APPENDIX  A. 


FORMULATION  OF  ACOUSTIC  REFLECTION  COEFFICIENT  AND  IMPULSE  RESPONSE 


Introduct i on 


A.  I 

This  Appendix  summarizes  part  of  the  analysis  needed  to  support  the  COAT 
project.  Although  the  results  were  derived  for  COAT,  they  are  quite 
general  and  applicable  to  many  other  problens.  The  general  expression  for 
the  acoustic  plane  wave  reflection  coefficient  of  a  plane  layer  is  derived 
(using  a  wave  number  Hatching  approach)  and  applied  to  a  viscoelastic 
layer.  The  results  of  this  derivation  have  been  inplemeoted  in  the  FORTRAN 
code  program  COAT  which  calculates  the  total  reflection  coefficient  as  a 
function  of  layer  thickness,  incidence  angle,  frequency,  and  material  pro¬ 
perties  (at  a  given  temperature). 

The  layer  (reflected)  impulse  response  is  then  developed  in  general  form 
and  application  to  real  problems  of  viscoelastic  layers  is  discussed.  The 
technique  has  been  implemented  in  the  FORTRAN  program  IMPULSE  which  uses 
the  reflection  coefficient  lata  determined  oy  COAT  to  calculate  the  impulse 
response  for  a  viscoelastic  layer.  The  results  from  IMPULSE  are  approxima¬ 
tions  of  the  actual  response,  and  the  effects  of  the  finite  EFT  approxima¬ 
tion  used  must  he  analyzed. 

The  results  of  applying  the  technique  on  an  elastic  layer  are  discussed  in 
an  example.  The  COAT  program  provides  the  reflection  coefficient  data 
needed  as  input  into  the  IMPULSE  program.  The  example  shows  excellent 
agreement  with  theoretical  predictions  for  reflected  pulse  spacing  and  good 
agreement  for  the  reflected  pulse  amplitudes  accounting  for  corrections 
associated  with  the  finite  FFT  technique. 
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Parameter  List 


A.  2 


Fluid: 

density 
sound  speed 


Viscoelastic  Layer: 

density  p2 
sound  speed  C2 
thickness  H 
dynamic  shear  modulus  p 
Poisson's  ratio  v 
shear  damping  factor  5 


complex  Young's  modulus 


Bigid  (Base)  Material  : 

density 

sound  speed  C3 


Acoustic  Plane  Wave: 

wave  number  k=;o/ c 

radian  frequency  w 

angle  relative  to  boundary  9 
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A. 3  CW  Reflection  Coefficient 

A  general  expression  is  derived  for  the  plane  wave  (pressure)  reflection 
coefficient  of  a  plane  layer  separating  two  other  media.  This  general 
reflection  coefficient  is  expressed  is  a  function  of  frequency  j  incidence 
angle  q  ,  and  material  properties  {■>.,€■).  This  -general  result  is  then 
applied  to  the  case  of  a  plane  / i scool astic  layer  on  a  perfectly  reflecting 
base  material  in  water.  The  water  is  designated  as  median  1,  the  visco¬ 
elastic  layer  as  medium  2,  and  the  base  oaten al  as  medium  3  (figure  A-l). 


A. 3.  I  Derivation  of  the  General  Expression: 


3y  examiniig  the  reflection  and  refraction  of  the  plane  wave  at  the  two 
boundaries  (13  and  23)  between  the  three  media  (Figure  A-l),  the  total 
plane  wave  reflection  coefficient  for  the  layer  is  found  to  be 


R  =  Rl2  +  F12^23T21  e  *  +  T12R23R2l'*?3T21  ~ 


where  =  boundary  reflection  coefficient  for  pressure  wave  traveling 

fron  median  i  towards  medium  j 

T^j  =  boundary  transmission  coefficient  for  pressure  wave  traveling 
from  medium  i  to  median  j 

0  =  Hk-,  cos  ^2  represents  the  phase  shift  in  the  layer  of  thick¬ 

ness  H  (where  k^  cos  09  is  the  vertical  component  of  wave- 
number  in  the  layer). 

Notice  that  the  expression  for  R  is  in  the  form  of  a  geometric  series  where 

l  r11  =  ( 1-r) 
n  =  0 
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so  that  the  total  reflection  coefficient  becomes 


R  =  Rl2  +  T12  R23  e-12^  [  l  R„  e-,2^)n]  T, 


21  23 


a  R12  + 


T  R  T  e"i2d 
12  23  21  9 


1  '  R21R239 


Since  T.^T^.  =  1-R2  and  Rp  =  _R21’  tf,e  e*Press’on  therefore 


9  +•  P  p 

12  23 


1  r  R12R23  e 


This  is  the  general  expression  for  the  plane  wave  (pressure)  reflection  coef¬ 
ficient  of  a  plane  layer  separating  two  other  media.  R  is  a  function  of  inci¬ 
dence  angle  0^,  frequency  w,  and  naterial  properties  (ppc^  »p? »c? ,p3 ,c3) 
through  the  expressions  for  phase  shift  <*>  and  boundary  reflection  coefficients 
R32  and  R23  given 


♦  =  Vn12  '  Sln2  *1 


m.  .  cos  0i  -  \/n?.  -  sin2e. 

^  —  '  J  ^  *  J  ' 

1J  m..  cos  9.  +  Vn  —  -  sin2  0. 
ij  i  v  ’J  i 
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I 


where  H 


=  layer  thickness 


in:  :  =  p./Pi  density  rati o 

J  J  * 

n^j  =  c-j/cj  index  of  refraction 

i  =  1,2 

j  =  2,3. 


A. 3. 2  Application  To  A  Viscoelastic  Layer 

In  the  case  of  a  viscoelastic  layer  on  a  rigid  base  material  where 
p?c^/p3C3<<l ,  ^23-+l  (hard  boundary  with  no  phase  inversion  on  perfect 
reflection).  For  viscoelastic  materials,  the  compressi onal  wave  (sound) 
speed  c?  is  a  function  of  'naterial  properties  (which  are  dependent  on 
frequency  <*>)  such  that 

r  =  r  (1-v)  E  -1/2 

2  : (  l  +  v) ( l-2v)  p?J 

(notice  that  since  c?  =  C2(w),  the  material  is  dispersive).  Therefore,  for 
a  viscoelastic  layer  separating  water  and  a  perfectly  reflecting,  rigid 
base  material,  the  plane  wave  reflection  coefficient  becomes 


n  ~i 

p  +•  e  * 

3  12 

1  *  ri2  e"24 

where  .{>  =  H  —  D? 

l 

phase  shift  term 

D1  "  n2 

12  f 

boundary  reflection  coefficient 
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c 


2 


r  ( 1  -  v»)  E 
'  (l  +  v) ( l-2v)p2' 


compressi onal  wave  speed 


E  »  Er  (Hi5) 


complex  Young's  modulus 


Er  =  2u  (1+v) 


real  part  of  Young's  modulus 


v  =  constant 


Poisson's  ratio 


P  =  p(w) 


dynamic  shear  modulus  (real) 


5  =  5(o) 


loss  factor 


(notice  that  3  is  complex  and  accounts  for  absorption  and  attenuation  due 
to  the  viscoelastic  layer). 


A. 4  Inpul se  Response 


A  general  expression  of  the  impulse  (reflected)  response  of  a  plane  layer 
for  an  incident  plane  acoustic  impulse  is  derived  using  the  CW  reflection 
coefficient  of  the  previous  section.  A  method  for  finding  this  (reflected) 
response  for  viscoelastic  layers  in  real  conditions  is  then  discussed. 


A. 4.1  General  Impulse  Response  Integral 

The  incident  pressure  field  of  a  harmonic  plane  wave  of  unit  amplitude  can 
be  written 

ik.(xsine  -zcosQ  )  -  iwt 

P.  =  e  1  1  1 

me 


and  the  corresponding  reflected  field  is 


?refl 


ik.(xsino.  +  zcose.)  -i-ot 
Re1  1  1 


(Figure  A-2)  where  R  is  the  total  reflection  coefficient  (given  in  the  pre¬ 
vious  section).  The  impulse  response  (general)  is  found  by  integrating  the 
reflected  field  over  all  frequencies  so  that 

.  +<=  ik,(xsine,  +  zcose,)  -i-ot. 

I  -  -k  I  «  e  1  1  1 

-  00 

and  since  =  oj/c^,  the  response  integral  can  be  expressed  as  an  inverse 
Fourier  transform  of  the  form 

I  (x)  =  -L  J  R(uj)  eiuxdu, 

—  CO 

where  R(w)  =  CW  reflection  coefficient  (previous  section) 
t  -  (xsino.  f  zcosQ,  -  c.t). 

C  i  1  L  L 


A. 4. 2  Solution  For  A  Viscoelastic  Layer 


To  solve  the  response  integral  I(-c)  for  actual  material  in  real  conditions, 
the  FFT  method  appears  to  be  the  most  efficient  since  an  analytical  solu¬ 
tion  is  not  probable  due  to  the  complex  functional  dependence  of  R  on  u 
(described  in  the  previous  section).  The  frequency  dependence  of  the  layer 
material  properties  n(/j)  and  6(m)  are  critical  in  determining  whether  or 
not  the  integral  I(x)  is  integrable  in  its  present  form  and  what  minimum 
and  maximum  frequency  limits  are  appropriate  for  the  FFT  to  obtain  an  ade¬ 
quate  approximation  of  the  integral.  Also,  since  the  integration  is  taken 
from  negative  to  positive  frequencies,  and  the  material  properties  are 
defined  only  for  the  positive  frequencies,  values  of  R(w)  must  be  extrapo¬ 
lated  for  negative  frequencies. 

In  order  to  find  the  values  of  R(<a)  for  negative  frequencies,  the  observa¬ 
tion  that  (in  reality)  the  impulse  response  is  real  only  (not  complex)  is 
used  to  determine  the  .-ven/odd  characteri sti cs  of  R (m).  The  integrand  of 
I ( t)  is  a  complex  expression 

Re1jr  =  Re  { Re 1  0T }  +  i  Im  (Re1  n]  . 

For  the  integral  I ( x)  to  be  real  only,  Im  {RelwTl  must  be  an  odd  function 
(whose  integral  from  -°°  to  <-=>  is  ^ero).  Since  Re(e1ux}  is  even  and 
IM { e1 }  is  odd,  then  Re{R|  must  be  an  even  function  and  Im(R}  must  be  an 
odd  function.  Therefore,  the  value  of  R(to)  for  negative  frequencies  can  be 
extrapolated  from  the  results  of  R(w)  for  positive  frequencies  by  applying 
the  knowledge  that  the  real  part  of  R(«)  is  an  even  function  and  the  imagi¬ 
nary  part  is  an  odd  function.  This  can  be  expressed  as 

Re { R ( - o) }  -  Re{R(m) }  even 

Im{R(-<o) }  =  -  Im{R((j) }  odd 

or  R(-m)  =  R*( w) 
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where  the  *  superscript  denotes  the  complex  conjugate. 

The  minimum  and  maximum  frequency  limits  for  the  FFT  are  determined  by  the 
material  properties  data  available.  Extrapolation  of  the  material  proper¬ 
ties  data  beyond  the  frequencies  for  which  they  are  available  (when  greater 
accuracy  in  the  impulse  calculation  is  needed)  can  be  based  on  experimental 
observations  and  a  physical  understanding  of  the  important  processes.  The 
effect  of  approximating  the  impulse  response  integral  with  a  finite  FFT 
must  still  be  determined.  The  correction  coefficient  for  the  reflected 
response  to  a  unit  amplitude  impulse  is  dependent  on  the  frequency  limits 
used  in  the  finite  FFT  approximation. 


A. 5  Application  of  Technique 


The  technique  of  calculating  the  reflection  coefficient  (from  material  pro¬ 
perties)  and  the  reflected  impulse  response  of  a  plane  layer  has  been 
applied  to  the  case  of  an  elastic  layer  for  which  the  exact  theoretical 
solution  is  tractable.  In  this  way,  the  technique,  discussed  in  the  pre¬ 
vious  sections,  can  be  compared  with  the  exact  solution  tc  find  the  degree 
to  which  this  finite  FFT  technique  approximates  the  expected  reflected 
response. 

The  programs  COAT  and  IMPULSE  have  been  tested  for  a  10  cm  thick  elastic 

layer  with  E  =  2.83  x  10^  dynes/cm^  and  v  =  0.44  for  a  unit  impulse  at  nor¬ 

mal  incidence.  The  density  of  the  layer  was  taken  to  be  twice  that  of  the 
first  medium  (p^/p^  =  2).  The  results  (calculated  by  COAT)  for  the  real  and 

imaginary  parts  of  1  (the  reflection  coefficient  of  the  layer)  as  a  func¬ 

tion  of  frequency  (up  to  104  Hz)  are  plotted  in  Figure  A-3.  The  reflected 
response  to  a  unit  impulse  (calculated  by  IMPULSE)  is  plotted  in  Figure 
A-4.  The  F FT  size  was  2043  (for  a  full  frequency  range  of  10^  Hz)  with 
cosine  tapering  of  10%  on  each  end. 

The  exact  theoretical  prediction  for  this  normal  incidence  case  is  given  in 
Figure  A-5.  The  theoretical  prediction  gives  a  reflected  pulse  spacing  of 
about  1  msec  and  reflected  pulse  amplitudes  (ratios  of  reflected  and  inci¬ 
dent  amplitude)  of  -0.561,  0.685,  0.384,  0.216,  0.121,  0.068,  0.038,  etc., 
respectively.  Comparing  the  expected  (theoretical)  result  with  the  IMPULSE 
calculation  (Figure  A-4),  we  see  that  the  (IMPULSE)  technique  correctly 
predicts  the  pulse  spacing  and  closely  approximates  the  amplitudes  of  the 
odd  pulses,  but  the  even  pulses  (starting  with  the  first  positive  puise) 
show  much  less  agreement  in  amplitude.  Increasing  the  size  of  the  FFT  (the 
frequency  resolution  of  the  data)  did  not  significantly  affect  this  discre¬ 
pancy  between  odd  and  even  pulse  amplitude  approximations ,  but  increasing 
the  frequency  range  (and  limit  of  the  FFT)  up  to  105  Hz  eliminated  this 
discrepancy  at  the  expense  of  increasing  the  error  for  the  odd  pulses 
(Figure  A-6). 
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The  results  of  this  example  indicate  the  effects  of  approximating  the 
impulse  response  integral  with  the  finite  FFT  technique  and  possible 
methods  of  correcting  for  the  discrepancy  (i.e.,  increasing  the  frequency 
limits  to  reduce  the  odd/even  pulse  error).  The  COAT/IMPULSE  technique  's 
shown  to  be  a  good  predictor  of  the  expected  (actual)  reflected  impulse 
response  given  the  geometry  and  material  properties  of  the  layer  of 
interest. 
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Figure  A-l„  Plane  layer  (plane  wave)  reflection. 
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Figure  A--2.  Schematic  of  plane  layer  total  reflected  layer. 
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